Field induced d x 2_ y 2 + id xy state in d-density-wave metals 
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We argue that the d xy component of the order parameter can be generated to form the d x 2_ y 2 + 
idaiy-density wave state by the external magnetic field. The driving force for this transition is the 
coupling of the magnetic field with the orbital magnetism. The fully gapped particle spectrum and 
the magnetically active collective mode of the condensate are discussed as a possible signature of 
the d + id' density wave state. 
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Recently, Chakravarty et al. Q proposed to model the 
pseudogap (PG) of underdoped high-T c cuprates in the 
formation of a new order — d-density wave (DDW) state, 
which breaks the parity and time-reversal symmetry, and 
the invariance of translation by one lattice constant and 
rotation by ir/2. The scenario is supported by the experi- 
ments showing that the PG and the superconducting gap 
coexist distinctly below T c . Although the debate about 
the mechanism for the PG is far from settled, it is sugges- 
tive that the typical charactistics of the PG as observed in 
many experiments, including photoemission tunnel- 
ing |3j , muon spin relaxation Q , can be explained by the 
DDW model. More recently, the detection of the DDW 
ordering using impurity resonance has been proposed || . 

We are going to argue here that, in addition to the 
dominant d x 2_ y 2 (d) component of the DDW order pa- 
rameter, a subdominant d xy (d') component can be gen- 
erated by the magnetic field. We find that i) the pres- 
ence of the additional gap component d' will lead to the 
fully gapped quasiparticle spectrum; ii) the inherently 
two component order parameter d + id' of the ddw phase 
in the field will exhibit collective mode that can be ex- 
cited in an out-of-plane ac magnetic field. Both of these 
features can be used to distinguish between ddw scenario 
and alternatives. We discuss experimental consequences 
below. The existence of the field induced d + id' super- 
conducting state in high-Tc cuprates has recently been 
intensively studied. 

The existence of the field induced d + id' supercon- 
ducting state in high-Tc cuprates has recently been in- 
tensively studied. It was first proposed by Laughlin 
and Ramakrishnan JtJ, to explain the kink behavior ob- 
served in the thermal transport experiment [||, that the 
magnetic field can drive the original d-wave ordering into 
the d + id' state. In the original proposal, the generation 
of the d! component was thought to be through the first- 
order bulk phase transition, but the important role of the 
vortices for the anomaly was also pointed out by other 
authors || [h], . Later on it was argued Q that even 
when the magnetic field is near the upper critical field, 
the d-wave state can be "distorted" by the external field, 
leading to a bulk d + id' -wave state with the intrinsic or- 
bital moment. There is a significant distinction between 



the d-wave superconducting (DSC) state and the DDW 
state: The magnetic induction in the superconducting 
state is inhomogeneous. In the DSC state, the electro- 
magnetic f/(l) gauge invariance is broken so in an exter- 
nal magnetic field, the system is in either the Meissner 
state (type-I superconductors), where the magnetic flux 
is expelled from its bulk region, or nucleating an array 
of vortices (type-II superconductors). The DDW state 
does not break the gauge invariance and has no analog of 
the Messner effect or the Abrikosov vortices. It is then 
expected that magnetic field inside a sample in the DDW 
state is homogeneous. In view of this distinction, we ar- 
gue that there should be an instability of the d-density 
wave state into the d + id'-density wave state even in the 
presence of a weak magnetic magnetic field. 



As in the DSC state, there also exists low-lying quasi- 
particle states around the nodes of the <i-density wave 
energy gap. At these nodes, the low-lying quasiparticles 
have vanishingly small energy gap and they are most sen- 
sitive to an external electromagnetic perturbation, which 
opens the possibility of the generation of the second com- 
ponent of the order parameter, orthogonal to the initial 
d-density wave state. Second component of the order 
parameter can be of s-wave or d'-wave orbital symme- 
try. Although, the s-wave component may be induced 
due to the scattering at surfaces or interfaces, the d'- 
wave order parameter is likely to be generated in the 
bulk sample when an external magnetic field is applied. 
The physical origin of this instability is the bulk orbital 
magnetic moment (M z ) in the d + id' state. When an 
external magnetic field H is applied perpendicular to the 
plane of the two-dimensional (2D) system under consid- 
eration (namely, H || z), the resulting coupling of the 
magnetic induction B with the orbital magnetic moment, 
— {M Z )B, lowers the system free energy. As mentioned 
above, since in the DDW state there is no screening ef- 
fect on the magnetic field, the magnetic induction B is 
homogeneous throughout the system and is close to the 
external magnetic field H . In the absence of the mag- 
netic field, the pure d-density wave state can be regarded 
as the equal admixture of the orbital angular moment 
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L z = ±2 pairs: 

W (G) = iWocos(26) = ^[exp(2i9) + exp(-2i9)] . 

(1) 

Here we have made an approximation to the order pa- 
rameter Wo(k) cx (c^q .Ck, CT ) oc Wo(cosk x a — cosk y a) 
by confining the wave vector k near the Fermi surface and 
introduced as the 2D azimuthal angle of the Fermi mo- 
mentum, where c ki<T annihilates an electron of spin a at 
k, Wq is the magnitude of the pure d-wave component. In 
the presence of an external magnetic field, the L z = ±2 
orbital wave functions becomes unequal and the coeffi- 
cients for them are shifted linearly with the magnetic 
field H: 

W {&) -> ^[(l + 77S)exp(2ie) + (l-r ? B)exp(-2ie)] 
= i[W (Q)+iBWi(e)] , (2) 

where B = H and Wi » 77sin(26). Unlike the DSC, the 
pure (i-density wave order parameter is imaginary while 
the field generated d'-wave component is real. We will 
still call the resulting DDW state as having the d + id' 
symmetry because the relative phase between the two 
components is still ir/2 in the equilibrium. 

To support the above physical intuitive, we provide 
a microscopic analysis in the following. The d-density 
wave ordering comes from the formation of the particle- 
hole pairs caused by the electron interactions. The col- 
lective motion of these pairs can be represented by the 
center-of-mass coordinates, while the relative motion by 
the relative coordinate. The structure of the pair wave- 
function is determined by the relative motion of two 
paired particles. In the tight-binding model, both the 
DDW term and the coupling term between the magnetic 
field with the electron orbital angular momentum (i.e., 
the orbital Zeeman coupling) are associated with the bi- 
linear operator c\ a Cj^ a . Therefore, in the presence of 
magnetic field, both of these two terms should be as- 
signed, together with the kinetic energy term, a phase 
factor exp(^ J^' A • dl) to ensure the gauge invariance 
of the system Hamiltonian, where <i>o = hc/e and A is 
the vector potential. This fact indicates that the DDW 
order parameter is not influenced by the gauge phase fac- 
tor, and allows us to treat the orbital Zeeman coupling 
effect on the DDW order parameter independently. Ex- 
plicitly, the coupling of the magnetic field (parallel to z 
direction) with the orbital anugular momentum can be 

written as: H B = -*5Mb-B c i,a^( T i x $)zCi+s,<r- 

Here S = (±a, 0) or (0, ±a) is the unit vector along the 
x or y axis, a is the lattice constant, r< is the lattice vec- 
tors, /Lis is the Bohr magneton, g = 2m e a 2 t / ti 2 with t the 
nearest-neighbor hopping integral [the energy unit used 
thereafter] , m e the effective mass of electrons and H the 
Planck's constant. In the momentum space, the system 



Hamiltonian can be written as: 

H = 5ZCkC k)CT c kiCT + ^[W k c k)CT c k+Q , (T + H.c.] 

-igH B B } j c k a [sin kg x <9 kQ ] 2 c kj(T . (3) 

k,cr 

Here £ k = —2 [cos k x a + cos k y a] — ji with ji the chemical 
potential is the single particle energy measured relative 
to the Fermi energy. The DDW order parameter is given 
by Wu — zWoyo(k) + Wi<^i(k), where y>o(k) = cos k x a — 
cosk y a and yi(k) = sin fc^a sin fc^a. The amplitude of 
the d- and d'-wave compoponents Wo a are determined 
self-consistently: 

^ = ^E( c U,^)^o(k), (4) 

k 

and 

2U 

Wi = -^$>k +Q , CT c k , CT )^i(k) , (5) 

k 

where Vb,i are, respectively, the d- and <i'-channel inter- 
action, N is the number of 2D lattice sites. That the 
d-wave component is imaginary is due to the equiva- 
lence of Q = (tt,tt) and — Q enforced by the underly- 
ing band structure. The notation [sinka x dk a ]z rep- 
resents sin(k x a)dk y a — sm(k y a)dk^a- We define e k = 
— 2[cos k x a-\-cos k y a] so that £ k = e k — ji. For Q = (7r,7r), 
we have following symmetry properties: e k+ Q = — e^, 
<^ (k + Q) = — <y3o(k), and y>i(k + Q) = ^i(k). In view 
of the fact that the DDW state breaks the translational 
symmetry with lattice constant but conserves that by 
\/2a along the diagonals of the square lattice, it is con- 
venient to halve the Brillouin zone, by introducing two 
kinds of electron operators c k(T and c k +Q, CT . The pair- 
ing of the particles and holes must cause correlations in 
their relative motions. According to the structure of the 
Hamitonian and the self-consistency conditions for the 
DDW order parameter, we can introduce the following 
Green's functions to describe the correlation: 

0n(k,k';r) - -(^[^(^4,^(0)]), (6a) 

S 12 (k,k';r) = -(r T [ Ck+Q!(T (r)4, ^(0)]) , (6b) 

&i(k,k';r) = -(T T [c Ka (r)cl +Q jO)}) , (6c) 

S 22 (k,k';r) = -(r r [ Ck+Q , ff (r) Ck , +Q!(T (0)]) , (6d) 

where the factor T T is a r-ordering operator as usual, 
c\l.<j{t) — e Hj c\ L a e~ Hr is the operator in the Heisenberg 
representation. Given the Hamiltonian Eq. (|^), with aid 
of the equation of motion for the field operator c k ^(r) 
and c k (t), and by performing a Fourier transform with 
respect to r, 

£(k,k';r) = fc S T^£(k,k'^„)e- M " T (7) 
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with u> n = (2n+ l jirksT, we establish a closed set of the 
self-consistent equations for G(k, k'; iu) n ), e.g.: 

<5kk' = (ico n - Ck+Q + i#Ai£!-B[sin(k + Q)a x d( k+Q ) a ] z ) 
x£ 22 (k, k'-iu n ) - 2{W£ + VF k+Q )£ 2 i(k, k'; iu n ) , 

(8a) 

and 

= {iuj n - £ k + i^s-B[sinka x <9 k a]z) <?2i(k, k'; iuj n ) 
-2(W k + W k * +Q )fe(k, k'; iw n ) . (8b) 

To the approximation up to the first order 
in the orbital-magnetic field coupling, we obtain 
£/2i(k, k'; iuj n ) = ^(k, k'; iw n ) + 50 2 i(k, k'; zw„), where 

ft(k,^- w ^r^ , (9, 

-igfi B B(iuj n - Ck+Q) 

D(k; iw„) 
x[sinka x <9 kQ ]£ 2 i (k, k'; iw„) (10) 



and 



<50 2 i(k,k';iw n ) 



where _D(k;iw„) = (iu> n ~~ E^i)(iuj n — -Ek.2) with 



#k,i(2) = ±yJ<i + \Ww + W£ +Q \ 2 - pl. We take the 

ansatz that Vq is bigger than V\ jl3| such that in the ab- 
sence of the magnetic field, the (i-wave ordering is pure 
and no secondary phase transition for the appearance of 
the dl ordering occurs. Therefore, the DDW gap appear- 
ing in the Q° is, Wk = ?W / o < /?o(k). The momentum de- 
pendence of <£>o(k) leads to [sinka x 9ka]<^o(k) = 2ipx(k). 
Substitution of Eq. @ into Eq. @ yields: 



Wi 



4Vi 

W 

r]BW 



Re[50 2 i(k,k;r = O)]^(k) 



kerbz 



(11) 



where 



16g^BVifc B r 
N 

16 gM BiV(0)Vi 
E F 



E E 



kerbz 



D 2 (k;iu) n 



(12) 



where N(0) is the density of states at the Fermi energy 
Ep. By taking the Fermi wave length of a few lattice 
constant a (~ 4A) and AT(0)|V1| - 0.3, it is estimated 
|Wi/Wo| « 10~ 2 at fl = 10T, which makes the ampli- 
tude of the induced component | W\ | to be on the order 
of a few Kelvin. Eq. (|ll|) shows microscopically that the 
magnetic field can drive the pure d-wave state into the 
d + id' state. This analysis also indicates that the d + is- 
density wave state could not be generated by the applied 



magnetic field, which is significantly different from the 
situations at the sample surfaces or interfaces. 

Equation ( pd| ) suggests that the phenomenological 
Ginzburg-Landau (GL) free energy functional must con- 
tain the linear coupling between the original d-density 
wave order parameter and the field-induced d'-density 
wave order parameter, / m ^. = i^(iWo)W*B + c.c. Con- 
sequenetly, we can write down the system GL functional 
of the form: 



d\\^{T-T« c )\W,{v)\ 2 



0o 



W (r)\ 4 



-^|V^o(r))| 2 + ^|VWMr)| 2 + y |Wi(r) 



+/int(r)l 



(13) 



where the first two terms describe the instability of the 
pure d-density wave state, with Xj? being the transition 
temperature in the absence of magnetic field. The last 
two terms represent the energy shift of the d-wave state 
as a result of the field-induced d'-wave order parameter, 
where ct\ is positive. Notice that, unlike the supercon- 
ducting order parameter, the gradient operator on the 
DDW order parameter is not shifted by the vector po- 
tential because the DDW pairs do not carry charge. It 
follows from Eq. ( p^[ ) that, as far as the d'-wave com- 
ponent is concerned, the coupling to the magnetic field 
term [i.e., / m ^] is linear, while the stiffness term [i.e., the 
second last term in Eq. (O)] is quadratic. Therefore, at 
least at the weak field so that W\ is small, the linear term 
is dominant. Therefore, the system gains energy by hav- 
ing a nonzero equilibrium value of W± . By treating Wi 
and W* as independent variables, the GL functional T 
is minimized by enforcing = = 0, which leads 
to 



riB 

Wi = —W 

Oil 



(14) 



Upon substituting the above result into Eq. (|13J) , we find 
the energy gained by the system with the induced d'- 
density wave component: 8 J- = — J d 2 rr) 2 \ Wq\ 2 B 2 /2a\. 
Therefore, the transition temperature which is now field 
dependent, and is renormalized by the magnetic field as: 
T C (B) = T° + ST C (B), where ST C (B) = r 1 2 B 2 /2a a 1 . 
It then follows that coupling of the magnetic field with 
the orbital angular momentum shifts the instability of 
the d-density wave ordering to the high temperature as 
schematically shown in the (T, B) phase diagram. Here 
we note that, since the particle-hole pairing takes place 
with the equal spin, the coupling between the magnetic 
field and the electron spin (i.e., the spin Zeeman cou- 
pling) will not depress the induction of d' component in 
the DDW metal. We do not address here the effect of 
very strong field that could ultimately supress the DDW 
state. 

Up to now our analysis of the induction of the sec- 
ondary dl component has been focused on the equilib- 
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FIG. 1: The schematic drawing of the {B,T) phase diagram 
of the d-density wave metal. The pure d x 2_ y 2 density wave 
state is marked by the thick solid line lying on the y axis. 
Quite different from the case of d-wave superconductors, the 



d. 



„2-density wave state has the instability into the d x 2_ y 2 + 



idxy-w&ve state in the presence of an arbitrary magnetic field 
coupled to the electron orbital angular momentum. 



rium solution. If we assume that this secondary d! or- 
der parameter has been created, we can write in general 
iW Q = |W |e^° and Wi = |Wi|e*, and study the dy- 
namics of the relative phase <p — 4>i — 4>o, which is gov- 
erned by El: 



d 2 ^ 
8t 2 



where p- 1 « N{0). With Eq. (03h, we find 



a 2 
dt 2 



-£ = -p- 1 r ] B\W \\Wi\cosc 



s 2 V 2 c 



(15) 



(16) 



which leads to the clapping mode with dispersion 
uP{B,k) = uj 2 {B) + s 2 k 2 with u)q(B) = r]B 2 \W a \ 2 /p and 
s 2 = \W \ 2 (K + r] 2 B 2 Ki) / 4p. This mode represents the 
oscillation of the relative phase between the d and d! com- 
ponents of the DDW order parameter, and is tunable by 
the magnetic field. 

We thus have proved that (a) the applied magnetic 
field can generate the d xy order parameter in the d- 
density wave metal, whose amplitude is linearly propor- 
tional to the field strength, (b) the transition into the 
d + id'-density wave state occurs at a higher transition 
temperature, and (c) there exists a new clapping mode 
corresponding to the oscillation of the relative phase be- 
tween the two components. We now turn to the exper- 
imental observation of the existence of the induced d' 
com ponent. Angle rcaolvcd photocmiaaion spectroscopy 



FIG. 2: Energy distribution curves at (-|, for various values 
of the d xy -w&ve component W\ = (solid line), 0.01 (dashed 
line), and 0.02 (dotted line). The other parameter values: 
p = and T = 0.005. 



surface while d'-density wave gap reaches the maximum. 
Figure |2| displays the low-temperatures ARPES signal 
J(k, u>) = f(uj)A(k.,uj) at k = (|, |) for various values of 
the induced d'-wave component W\. Here is the Fermi 
distribution function f(u>) — l/[exp{uj /UbT) + 1] and the 
spectral function j4(k, u>) = 2[6(uj — -Ek.i) + 8(u> + E^)]- 
As is shown, the spectral peak is shifted with the mag- 
nitude of W\, which is in turn linearly proportional to 
the magnetic field. As another consequence, in the pres- 
ence of Wi, the quasiparticle spectrum is fully gapped. 
Therefore, if the magnetic field is applied perpendicularly 
to the 2D DDW metal, the electronic specific heat at low 
temperatures would be exponentially decaying. 

To conclude, we show that the c?-density wave state has 
an instability into the d + id' ordering in the presence of 
the magnetic field. The field induced d'-wave component 
is proportional to the field strength. The mechanism for 
the induction of the d' component is purely the coupling 
between the magnetic field and the orbital anugular mo- 
mentum. 
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